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Relative branching fractions of s — » h + h ~ decays (where h,h' = K or n) and direct CP asymmetry 
in B^ — > K + tt~ have been measured with 180 pb — 1 of data collected by the CDFII detector at the 
Tevatron collider. This includes the first BR measurement of a charmless B s — >PP decay (B B — > 
K+K-). 



1 Introduction 

The comparison of measurements performed 
on a number of such modes at e + e~ experi- 
ments with theoretical expectations has pro- 
vided a wealth of information. 

The CDF II experiment 1 has access to 
a large yield of these decays at the Teva- 
tron collider, thanks to a dedicated trigger 
on secondary decay vertices . In addition, it 
is also sensitive to the corresponding decays 
of B s mesons and beauty baryons, which pro- 
vide much additional valuable information 2,3 . 
In particular, combining measurements of B s 
modes with measurements on the B u /Bd sec- 
tor allows eliminating or constraining uncer- 
tain hadronic parameters. 

In this paper we report on the analysis 
of 2-body charmless decays of Bd, B s and A& 
in a sample of 179 ± 11 pb _1 collected with 
the CDF II detector between February 2002 
and September 2003. 

2 Sample Selection 

The data sample was selected in the 3-level 
CDF trigger system by a set of requirements 
dedicated to B — ► h + h ~ candidates. Two 
oppositely-charged tracks are required, with 
p T > 2.0 GeV/c, total p T > 5.5 GeV/c, im- 
pact parameters larger than 100 /jm, and az- 
imuthal opening angle between 20° and 135°. 
The B candidate is then required to have a 



transverse decay length larger than 200 /xm, 

invariant mass between 4.0 and 7.0 GeV/c 2 , 

and to point back to the primary vertex 

within 140 fim. 

In offline analysis, after reconfirming the 

trigger cuts with offline quantities, further 

cuts are imposed on transverse momenta and 

impact parameters of the two tracks, and on 

transverse decay length, impact parameter 

and isolation of the B candidate. Isolation is 

defined as — pr(-g)^ — where the sum runs 
pt(B)+2_ i PT 

on every other track within a cone of radius 
1 in rj-tp space around the B candidate flight 
direction. This cut is particularly effective 
in rejecting combinatoric background. The 
set of cuts was chosen to maximize the quan- 
tity 5/(5 + B) 1 / 2 , where S is the number of 
signal events expected from detailed simula- 
tion, and B is the background estimated by 
extrapolating the sidebands of the data. 

The resulting distribution of invariant 
mass, with an arbitrarily chosen pion mass 
assignment to both tracks, shows a clear sig- 
nal of 893 ±47 events and a = 38 ±2 MeV/c 2 , 
with a peak S/B in excess of 2 (Fig. 1). 

Detailed simulation predicts sizeable, 
closely spaced signals in this mass region 
from two Bd modes: Bd — » 7r + 7r~, Bd — > 
K + tt~ and two B s modes: B s — > K~ir + , 
B s — ► K + K~, forming a single unresolved 
bump. In addition, the — * pir~,pK~ 
modes might appear as a slight excess around 



5.5 GeV/c 2 . 

3 Measurement of individual 
modes 

The relative contributions to the signal of 
each Bd,s — > h + h ~ component, and the 
CP asymmetry of the self-tagging Bd — > 
K + ir~ mode were measured by means of 
an unbinned likelihood fit which combines 
kinematics and PID information on the two 
tracks. The Likelihood for the i th event is 
written as: 

Ci = b- C bck + (1 - 6) • C stg 
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Figure 1. Afn-Tr distribution of the B — > h + /t ~ can- 
didates after all selection cuts. The result of the fit 
is overlaid. 



The index sig (bck) labels the contribution 
of signal (background), and b is the back- 
ground fraction. The signal likelihood func- 
tion is C sl 3 = J^j fj ■ Cf n ■ Cf ID , where the 
index j runs over all possible B e i. s ~ ► h + h ~ 
modes and the parameters fj are their rela- 
tive fractions to be determined by the fit. 

£km j g - m p r i n ciple a function of the in- 
variant mass of the track pair, with a mode- 
dependent mass assignment to each outgo- 
ing particle. In a rigorous approach , one 
would then need to write every term in the 
likelihood (including the background part) as 
functions of the joint (correlated) distribu- 
tionof all four possible masses. 

In order to simplify the problem, only 
two variables are used: (1) the invariant mass 
computed with the pion mass assign- 
ment to both tracks and (2) a signed momen- 
tum imbalance, defined as a = (1 — — ) ■ gi, 
where p\ (P2) is the smaller (larger) of the 
track momenta, and q\ is the charge of the 
smaller momentum track. It can be shown 
that the candidate mass calculated with any 
possible mass assignment to the two tracks 
can be written to a good approximation as 
functions of just M T7r and a, which therefore 
provide a compact summary of all available 
kinematic information (Fig. 2). 



The kinematic term is then written as: 

c Un = 2 \ ' J . P J a ) 

J crv^ ° 

where the Mj(a) are very simple ana- 
lytical expressions obtained by series expan- 
sion; detailed detector simulation shows that 
they are very accurate within the kinematic 
range of interest (Fig. 2). Pj(ce) is distri- 
bution of a for each signal mode after the 
effect of reconstruction cuts, as determined 
from full simulation and parameterized by 
a 6 th order polynomial. The mass resolu- 
tion a is set by rescaling simulation predic- 
tions to match the observed widths of sev- 
eral other two-body decays (D° — > h + h ~ , 
J/tp — > an( l T ~~ * A 1 The uncer- 

tainty in this rescaling is accounted for in the 
final systematics. 

The Likelihood term related to particle 
identification information is 

Cf ID =pdff ID (\D 1 ,\D 2 ,a 1 ,a 2 ) 

with the observable ID carrying information 
from the specific energy release of the track 
in the drift chamber gas dE j 'dx meas , and is 

dE _ dE 

defined as: ID = £[ m ° as *2£ v ~" where 

dx cxp — K dx exp — "K 

dE j 'dxexp-K (dE / dxexp-n) is the expected 
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Figure 2. (Left) Af^x vs a for — > tt+.K' , com- 
pared to the first-order expression: = M^ + 

(2 - a)(ml -m\). (Right) Same for B d ^ -n~ K+ . 

value in the kaon (pion) hypothesis. The ex- 
pected value of ID is then by construction 
equal to for a pion, 1 for a kaon. This pa- 
rameterization allows using a single observ- 
able for both the pion and kaon terms. 

The function pdf PID 

is a non-factorizable distribution of two cor- 
related Gaussians, with track-dependent res- 
olutions o\ , (72 j and a fixed correlation a%i , 
determined from data. The distributions of 
resolutions have been checked to be equal for 
all modes 4 . The response of the drift cham- 
ber has been calibrated with a large sample 
of D° — > Kit decays from the same trigger, 
with their sign tagged by the presence of a 
D*^. The correlation between the measured 
dE/dx of two tracks in the same event, due 
to time-dependent fluctuations of the drift 
chamber gain, has been measured with the 
same sample. The average K/it separation 
was measured to be lAa in the kinematical 
range of interest. This moderate resolution 
is however sufficient to provide a statistical 
separation between K and ir which is 60% of 
what a perfect PID would provide. The un- 
certainty on the calibrations, and the effect 
of small unmodeled non-Gaussian tails have 
also been evaluated from data and included 
in the final systematic uncertainty. 

The background is described by C = 
£-bck ' ^6cl D ' with the kinematic term: 



where Pbck{ot) is the distribution of a for 
background events, obtained from sidebands 
of real data and parameterized as a 6 th de- 
gree polynomial. The Ci are free parameters 
in the fit. A few alternative parameteriza- 
tions of the background mass spectrum have 
been tested, and the corresponding parame- 
ter variations included in the systematics. 

The PID term for the background is simi- 
lar to the signal, and assumes that only pions 
and kaons are present, with each track having 
an independent probability f n to be a pion. 
The pion fractions f n are left free in the fit 
to vary independently in three mass regions 
(left, under, and right of the signal peak). 

The complete likelihood fit has been 
tested on MonteCarlo samples, and showed 
gaussian pulls on all variables with unit sigma 
and negligible biases. 

4 Efficiency corrections 

To extract relative branching fractions from 
the raw fit results it is necessary to correct for 
different efficiency of the selection on each de- 
cay. Most corrections are obtained from de- 
tailed detector simulation; exceptions are the 
efficiency of the B isolation cut (only affect- 
ing Bd/B s ratios), the difference in efficiency 
of the trigger track processor for kaons and 
pions due to their different energy release in 
the drift chamber gas, and the charge asym- 
metry. These three effects are not reliably 
simulated, and have been measured from real 
data. A systematic has been added for pos- 
sible differences between B d and B s momen- 
tum spectra 11 , that are treated as equal in 
the simulation. 

Isolation efficieny has been measured as 
a function of B transverse momentum from 
fully reconstructed samples of B s — ► J/V></>, 
B s D s n, B d -> J/i>K*, B d D-n+. 
The difference in trigger efficiency for pions 
and kaons has been measured from an unbi- 
ased sample of pions and kaons from charm 
decays. The momentum dependent charge 



asymmetry of CDF tracking has been mea- 
sured on large samples of inclusive tracks and 
on K s — > 7r + 7r~ decays 5 , and used to correct 
our measurement of direct Aqp. The size of 
this correction is about 1% and its relative 
uncertainty (±25 %) is included in the final 
systematics. 

The resulting total efficiency varies be- 
tween modes by less then 12%. 

The total systematic uncertainties are 
dominated by uncertainties in dE/dx cali- 
bration and efficiency of the isolation cut. 
Both uncertainties are due to statistical un- 
certainties in the calibrations samples, and 
are therefore expected to decrease with sam- 
ple size. 

5 Results 

The ratio of the Bd modes: 
BR(B d -> n ± n+) 



BR(Bd -> K+n*) 



0.24 ±0.06 ±0.05 



is in good agreement with the current world 
average 6 : 0.25 ± 0.025. The direct CP asym- 
metry is: 

N(lf d -> K-TT+) - N(B° d -> K+n-) 

N(B° d -> K-7T+) + N(B° d - K+n-) 

= -0.04 ±0.08 ±0.006 

which is compatible both with zero, and with 
the measurements performed by BaBar and 
Belle collaborations 8 . 
For the £?, sector: 



f s ■ BR(B S -> K+K+) 
f A ■ BR(B d - K±n+) 



0.50 ±0.08 ±0.07 



Since this is a CP eigenstate, it is possible for 
its lifetime to be different from the average B s 
lifetime measured from semileptonic modes. 
The above result is obtained under the as- 
sumption that T s = Td, the B s — > K + K~ 
mode is dominated by the short-lived com- 
ponent , and that Ar s /T s = 0.12 ±0.06 (the 
latter uncertainty is included in the quoted 



systematics). We also quote: 
fd-BR(B d ^TT ± n^ 



f s ■ BR(B S -» K±K*) 



= 0.48 ±0.12 ±0.07 



From the above numbers, the absolute 
BR can be obtained 10 : 

BR(B S -> ^K*) = 34.3 ± 5.5 ± 5.2 • 10~ 6 

This is almost twice the BR of Bd — ► K ± ir^ 
which differs only in the spectator quark, in 
agreement with recent calculations of flavor— 
SU(3) breaking 7,3 . This result will hopefully 
be a useful input for theoretical models of 
charmless decays. The large yield (w 230 
events in current sample) also offers an inter- 
esting opportunity for a measurement of Ar s 
without the need for an angular analysis. 

No evidence is found for the other B s de- 
cay for which a sizeable BR is expected, and 
a limit 12 is set: 

f a ■ BR(B S - K***) 



fd ■ BR(B d —* K+-K+) 
This translates to 10 : 



< 0.11 @ 90% C.L. 



BR(B S -> K^) < 7.5 • 10~ 6 @ 90% C.L. 

which is close to the lower end of current 
expectations. In addition to the above, rarer 
modes (BR: 10~ 8 to 10~ 7 ) dominated by 
penguin annihilation and exchange diagrams 
have been searched for by adding their ex- 
pected contributions to the Likelihood. In 
each case, the fit parameters changed by a 
negligible amount, and no evidence was found 
for those modes. The following limits 12 are 
set: 

BR{B d -» ^RT) 



or 



BR(B d -» K+nT) 

10. 



< 0.17 @ 90% C.L. 



BR(B d -» K^K*) < 3.1 • 10~ 6 @ 90% C.L. 

The current best limit in this mode is 9 0.6 • 
lO- 6 . 

BR(B S -> 7r ± 7 r : F) 



BR(B S -» K±K+) 



< 0.10 @ 90% C.L. 



This is derived under the assumption that 
both modes have the same average lifetime. 
This translates 10 to: 

BR(B S -» n ± nT) < 3.4 • 1(T 6 © 90% C.L. 

which is a substantial improvement over the 
previous best limit 9 of 1.7 ■ 1CT 3 . 

The same data have also been used in a 
separate search for the charmless decay A" — ► 
pir~ ,pK~ , by counting the number of events 
found in the search window 5.415 — 5.535 
GeV/c 2 . The choice of the window and 
the cuts have been optimized for maximum 
sensitivity 13 . No evidence for signal is found, 
and a Bayesian upper limit, based on uni- 
form prior, is obtained at 90% credibility: 
BR(A° b -> pn-fpR-) < 22 x 10~ 6 . 

This is an improvement over the previous 
limit: BR{h? b -» p-K~ ,pK~) < 50 x 10~ 6 at 
90%CL 14 . Predictions 15 lie in the range 2 — 
3 x 10~ 6 . Improved sensitivity is expected in 
the future with inclusion of PID information 
for the proton. 
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